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CN) ; ABSTRACT 

We present the results of fitting ASCA spectra of six narrow-line Seyfert 1 (NLS1) 
galaxies with the ionised reflection models of Ross & Fabian (1993). We find that five 
of the galaxies (TON S 180, PKS 0558-504, Ark 564, Mrk 335 and PG 1244+026) are 
well fit by the ionised disc model, and these are often better fits than the alternative 
models considered. The sixth galaxy, NGC 4051, has additional spectral complexity 
C\| | that cannot be well described by a simple ionised disc model or any of the other al- 

ternative models. The highest luminosity NLS1 considered, PKS 0558-504, does not 
have a well constrained ionisation parameter or reflection fraction. This is because it is 
difficult to distinguish between highly ionised, highly reflective discs and moderately 
ionised discs with low reflection fractions. The four galaxies with well constrained fit 
parameters are consistent with having inclination angles lying between 15 and 30 de- 
grees. Furthermore, these four sources are consistent with having a disc emissivity law 
that varies as r~ 25 . These last two properties are also typical of broad-line Seyfert 1 
| galaxies. We find little or no indication of a correlation between the reflection fraction 

. and the photon index of the underlying continuum. All six of the NLSls we considered 

show evidence for a broad Fe Ka line, but none of the line centroids are consistent with 
emission from highly ionised Fe. This is most likely due to the line being redshiftcd 
because of relativistic effects. We note that sources with larger ionisation parameters 
tend to have larger Fe Ka EWs. We interpret this as evidence that ionised Fe features 
are making their presence felt in the spectra. Since most of our sources have steep 
spectra, highly ionised features are predicted even by the new variable density reflec- 
tion models. The one source we analyse with a photon index less than two has the 
lowest ionisation parameter in our sample. We conclude that our result is the strongest 
evidence yet that NLS1 might have ionised accretion discs. This result gives further 
rS ' weight to the hypothesis that these objects contain rapidly accreting black holes. 
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1 INTRODUCTION 

One of the most basic difficulties in attempting to under- 
stand the physics of Active Galactic Nuclei (AGN), is to ac- 
count for the wide range of properties (i.e., broad lines, radio 
properties, etc.) that are observed over the whole menagerie 
of AGN classes. Although it is quite likely that we are ob- 
serving basically the same phenomenon through different 
observing angles (e.g., Antonucci 1993), not all AGN can be 
unified through this picture. Certain classes of AGN seem 
to be extreme examples of the phenomenon, and cannot fit 
into the simple unified scheme. It is hoped that such objects 
would be dominated by one physical property in their cen- 
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tral engines, and would therefore be simpler to model than 
a "typical" AGN. An example of this type of extreme AGN 
is the Narrow Line Seyfert 1 (NLS1) class of galaxies. 

NLS1 galaxies, like most AGN, are categorized by their 
optical emission properties. A NLS1 galaxy has the following 
characteristics: FWHM H/3 < 2000 km s" 1 , [O III]/H/3 < 3, 
and strong Fe II emission (Osterbrock & Pogge 1985). These 
properties mean that a NLS1 has narrower permitted lines 
(i.e., broad lines) then a typical Seyfert 1 galaxy. These 
galaxies also exhibit unusual X-ray properties (Brandt 1999; 
Leighly 1999a; Leighly 1999b; Vaughan et al. 1999b). In 
the soft X-ray band (usually taken to mean the ROSAT 
band: 0.1-2.4 keV) these galaxies show very steep spectra 
due to a significant soft excess, with a power-law photon in- 
dex (r; where the photon flux oc E~ T ) sometimes exceeding 
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3 (Boiler, Brandt & Fink 1996). In the hard X-ray band, 
the steep power-law flattens out to a typical F » 2.1-2.4 
(Brandt, Mathur & Elvis 1997; Leighly 1999b; Vaughan et 
al. 1999b) , which is significantly steeper than the spectra of 
broad-line Seyfert 1 objects. The soft excess is well fit by 
a model that is suggestive of thermal emission from a disc. 
NLS1 also exhibit more rapid and extreme variability in the 
X-ray bands than their broad-line counterparts (Turner et 
al. 1999; Leighly 1999a). 

The model that is currently the most successful in ex- 
plaining these properties is one which proposes that NLS1 
galaxies have smaller black hole masses than typical Seyfert 
galaxies (e.g., Boiler et al. 1996; Laor et al. 1997). This would 
explain the optical line properties, because the material in 
the broad line region would have smaller velocities and so 
their lines will have a smaller Doppler width. However, since 
these galaxies have about the same luminosity as broad-line 
Seyfert galaxies, then they must be emitting a higher frac- 
tion of their Eddington luminosity. This implies a higher ac- 
cretion rate and a hotter accretion disc, which would cause 
the thermal emission of the disc to shift into the soft X-ray 
band. The steeper high energy T could then be explained by 
Compton cooling of the hard X-ray emitting corona by the 
soft emission from the disc. The extreme variability would 
just be a result of the smaller black hole mass because the 
primary emission region would be smaller; however, it is dif- 
ficult to rule out the possibility of a beamed component. An 
analogy between NLS1 and Galactic black hole candidates 
emitting in their soft/high state has often been stated (e.g., 
Pounds, Done & Osborne 1995), although its usefulness has 
been debated. 

Hard X-ray observations of NLS1 can test this model. 
Models of Fe Ka line emission at various accretion rates by 
Matt, Fabian & Ross (1993) have shown that at modest ac- 
cretion rates, the accretion disc will have an ionised skin on 
top. This result assumes that the illuminating X-ray flux de- 
pends on the accretion rate, which is reasonable if the hard 
X-ray emitting corona does result from a process like mag- 
netic flaring (e.g., Galeev, Rosner & Vaiana 1979; Haardt, 
Maraschi & Ghisellini 1994; Svensson 1996) . Preliminary ev- 
idence for ionised accretion discs based on ASCA data has 
been found in two NLS1 (TON S 180: Comastri et al. (1998) 
and Turner, George & Nandra (1998); Ark 564: Vaughan et 
al. (1999a), but see Turner, George & Netzer (1999) for an 
alternative interpretation). However, these conclusions were 
made either by estimating the energy centroid of the ob- 
served Fe Ka line or Fe K edge and ignoring the rest of the 
continuum, or by using the pexriv ionised disc model of 
Magdziarz & Zdziarski (1995). Ionised reflection has observ- 
able effects over the whole continuum, so better constraints 
on the ionisation state can be made by fitting ionised disc 
models over the entire observed energy range. However, as 
shown by Ross, Fabian & Young (1999), the pexriv model 
is inaccurate as the disc becomes highly ionised. A further 
limitation of this model is that the Fe Ka line has to be 
added in separately. 

This paper presents the results of fitting the ionised disc 
models of Ross & Fabian (1993; hereafter RF) to ASCA ob- 
servations of 6 NLS1 galaxies. These models not only include 
the Fe Ka line, but also spectral features (emission lines and 
recombination continua) at lower energies. We take advan- 
tage of this information by fitting the models over a wide 



range of energy. The data selection and reduction procedure 
is described in Section bl We briefly outline the properties 
of our models in Sections, and then move on to report the 
results of our fits in Section ||. Finally, we discuss our results 
in Section p[ and summarise our conclusions in Section li. 



2 OBSERVATIONS AND DATA REDUCTION 

A summary of ASCA observations of the sample galaxies, as 
well as other relevant data, is shown in Table [y. The spectral 
data from the four detectors were reduced from the event 
files taken from the ASCA archive maintained by the ASCA 
Guest Observer Facility {ASCA GOF) at NASA/Goddard 
Space Flight Center, using FTOOLS 4.2 and the standard cal- 
ibration. Source photons were collected from a region with 
a radius of typically 4 arcmin for the Solid state Imaging 
Spectrometer (SIS; SO and SI) and 5 arcmin for the Gas 
Imaging Spectrometer (GIS; G2 and G3), while background 
data were taken from a source-free region on the detectors 
in the same observations. All the targets show significant 
X-ray variability during the observations. The ASCA light 
curves of our sample galaxies have been published elsewhere 
(Leighly 1999a). We present a spectral analysis of time- 
averaged data of individual sources integrated over each ob- 
servation. 

The SIS was operating with 2CCD chips switching 
between Faint and Bright modes for the observations of 
Mrk 335 and PG 1244+026. The Faint mode data of these 
observations have been converted to the Bright mode for- 
mat to add together with the original Bright mode data. For 
the others, the SIS was operating with 1CCD Faint mode 
throughout. These Faint mode data have been converted to 
a format called "Bright2" for analysis. 

Response matrices for the SIS were generated by SIS- 
RMG version 1.1. Version 4.0 matrices provided by the GIS 
team were used for the GIS. The effective areas of the source 
spectra were computed with ASCAARF version 2.73. For ob- 
servations carried out after 1994, the efficiency of the SIS, 
typically in the energy range below 1 keV, has been no- 
ticed to be decreasing continuously. Since the correction for 
this effect has not been implemented in ASCAARF at this 
stage, the effective area of the SIS in the low energy range is 
likely to be overestimated for some targets in our sample. We 
therefore restrict our spectral analysis to the energy range 
above 1 keV, except the data for Mrk 335 and NGC 4051 
which were observed when the SIS degradation was insignif- 
icant. 



3 IONISED DISC MODEL 

We employed the ionised disc models described in detail by 
RF and Ross et al. (1999), and so we will only give a brief 
outline here. A slab of gas consisting of solar Morrison & Mc- 
Cammon (1983) abundances, and with a constant hydrogen 
number density of uh = 10 15 cm -3 , is illuminated by an 
X-ray power-law with flux F x (defined over 0.01-100 keV) 
and photon index V. The computed reflection spectrum is 
then multiplied by a factor R, where R is the reflected frac- 
tion, and then added to the illuminating spectrum to find 
the model observed spectrum. 
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Table 1. Summary of ASCA observations and other relevant properties of the sample galaxies. Njj is the Galactic absorbing 
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Mrk 335 


0.025 


1993 Dec 9 


2CCD Faint/Bright 


25.1/20.4 


0.475/0.294 


1.50 


8.9 


3.04 


NGC 4051 


0.0024 


1994 Jun 7-9 


1CCD Faint 


81.7/71.1 


0.854/0.472 


1.31 


0.056 


2.84 


PG 1244+026 


0.048 


1996 Jul 1-3 


2CCD Faint/Bright 


49.6/38.8 


0.165/0.085 


1.93 


2.62 


3.26 


TON S 180 


0.062 


1996 Jul 10-11 


1CCD Faint 


54.2/49.2 


0.306/0.158 


1.50 


8.9 


3.04 


Ark 564 


0.024 


1996 Dec 23-24 


1CCD Faint 


54.8/50.8 


1.329/0.657 


6.40 


5.5 


3.47 


PKS 0558-504 


0.137 


1996 Sep 5-6 


1CCD Faint 


44.3/34.4 


0.660/0.359 


4.39 


124 


2.89 



Good exposure time of single detector of each type. 
2 Mean count rate in the 1—10 keV band obtained from the SO (for the SIS) and G2 (for the GIS) detectors. 
a data taken from Leighly (1999b). 



The important quantity in determining the structure of 
the reflected spectrum is the ionisation parameter, 



The larger £ is, the more ionised the gas is, and this will 
affect the strength and width of the features in the reflected 
spectrum like the Fe Ka line and the various absorption 
edges (Matt et al. 1993, 1996). 

We computed a three-dimensional grid of models with 
0.0 sC R < 2.0, 1.7 s£ T < 3.0, and 1.0 < logf < 6.0. The 
ionisation parameter was varied by changing only the inci- 
dent flux. We did not include the influence of soft radiation 
from the disc itself, but since we are mainly dealing with 
energies greater than 1 keV (Sections ^| & fl.l| ), this should 
not be an important effect. 

There are some important deficiencies in the RF mod- 
els. Foremost among them is the assumption of constant 
density in the reflecting slab. Recently, Nayakshin, Kazanas 
& Kallman (2000) presented ionised reflection models where 
hydrostatic equilibrium is solved along with the ionisation 
and radiation structure. They find that, due to a thermal 
instability, only a very thin layer (with Thomson optical 
depth, Tt Si 1) at the top of the slab is highly ionised, and 
the reflection spectrum is dominated by the cool, neutral 
material underneath. If this is the case, then it would be in- 
appropriate to fit data with the reflection spectra predicted 
by constant density models, such as the RF model used here, 
because they predict the presence of highly ionised features. 
However, all is not lost, because the models of Nayakshin et 
al. (2000) show that if the disc is radiation pressure domi- 
nated and the illuminating power-law has F <; 2 then ionised 
features can be seen in the reflection spectrum. In that case, 
our fits to these NLS1 sources might be reasonable because 
most of them have F > 2. Evidence for this interpretation 
is presented in Section H. 



4 RESULTS 

4.1 Spectral fitting 

Data from all four detectors were used in the spectral fit- 
ting. However, as mentioned in Section ^, the SIS detectors' 
efficiency below 1 keV has slowly been degrading since late 
1994. This is not corrected by the data reduction software 
and could lead to inaccuracies in determining continuum 
properties. Moreover, the sensitivity of the GIS detectors 



falls off sharply around 1 keV. To deal with these issues we 
adopted the following conservative approach. If an observa- 
tion was made after 1994 we determined if there was any 
discrepancy in fit residuals between the SIS and GIS detec- 
tors below 2 keV. If there was we used the GIS data down to 
1 keV and the SIS data only where the two agreed. If there 
was no visible discrepancy, both the GIS and SIS data was 
used down to 1 keV. If data was taken before late 1994 then 
the low energy SIS data may be used, but we made sure 
it was consistent with the GIS data. In all cases we used 
data up to 10 keV. The only exception to this was Ark 564, 
where we used simultaneous RXTE data kindly provided 
by S. Vaughan. This data was processed as in Vaughan et 
al. (1999a), except we used data from all 5 of the satellite's 
Proportional Counter Units. To avoid calibration and back- 
ground uncertainties we only fit these data between 3 and 
20 keV. Finally, we have noticed that the SO data of Ark 564 
are inconsistent in the energy range above 4 keV with the 
other four detectors which are in good agreement with each 
other. Since applying even more careful data selection than 
the normal practice could not eliminate the problem, we at- 
tribute it to some anomaly in the detector. We therefore use 
only the SO data below 4 keV in this particular data set. 

We used xspec vll.O for the spectral fitting. The grid 
of ionised disc models that was computed (Section ^) was 
converted into a tablefile that could be read in by xspec. 
In each fit, we let the relative normalisation between the 
SIS and the GIS vary, but we fixed the cold absorption at 
the Galactic value. Since we primarily fit data with energies 
greater than 1 keV, any excess absorption within a particular 
galaxy will not greatly affect our results. For each galaxy 
we fit its spectrum two different ways using the ionised disc 
models. First, we fit the data using just the model spectrum. 
Second, we included the effects of relativistic smearing by 
blurring the model with the Laor (1991) kernel for the Kerr 
metric. This blurring has four parameters: the inner radius of 
the emitting annulus (r m i„), the outer radius of the annulus, 
the exponent of the surface emissivity law (assumed to go 
as r a ), and the inclination of the disc along the line of sight. 
The outer radius was fixed to be at 1000 r g , where r g — 
GM/c 2 is the gravitational radius, (the actual value of this 
parameter will not affect the fit as the blurring effects come 
from the inner regions of the annulus), and the inclination 
angle was fixed at 30 degrees which is typical for Seyfert 1 
galaxies (Nandra et al. 1997a). We then fit the data with 
the models using four different combinations of r m i„ and a: 



© 2000 RAS, MNRAS 000, 0-| 



4 D. R. Ballantyne, K. Iwasawa & A. C. Fabian 



(r msn , a) = (10, -2.0; 10, -2.5; 6, -2.0; 6, -2.5), where r min 
is given in units of gravitational radii. Finally, using the best 
fitting model of the four considered, we fit the inclination 
angle. 

When the ionised disc models were blurred it was nec- 
essary for energies outside the response matrix of the de- 
tector to be used in order to avoid effects at the edges of 
the observed energy range. This was done by using the 'ex- 
tend' command in xspec to extend the energy range of the 
response matrix. As long as edge effects were avoided, the 
results of the fits do not depend on how far we extended the 
response matrix. 

The results of these fits are shown in Table ^ and in Fig- 
ures Plots of the ratio of the data to a simple absorbed 
power- law model can be seen elsewhere (Leighly 1999b). 
Note that all the spectra in Figs. |l]-^| are plotted in the 
observed frame, while all the fit parameters are quoted in 
the rest frame. The errors reported in Table [] and in all 
subsequent tables correspond to a 90 per cent confidence 
interval for one interesting parameter. 

To determine the Fe Ka line parameters of the galax- 
ies, we fit the spectra above 3 keV with a simple absorbed 
power-law model. To obtain a good estimate of the photon 
index, F, we ignored the data between 5 and 7 keV (observed 
frame) to avoid the influence of the Fe Ka line. Once we ob- 
tained the best-fitting value of F, we fixed it at that value 
and fit a Gaussian profile to the Fe Ka line. The energy 
centroid of the Gaussian was constrained to lie between 5.5 
and 7.5 keV, and its width (a) between and 2 keV. The 
rest frame fit parameters are shown in Table ^| Despite the 
differences in data reduction and analysis methodology, our 
values compare well with previously published results (e.g., 
Leighly 1999b; Vaughan et al. 1999b). 

As a check on the relative importance of our ionised 
disc fits we also fit the spectra with three alternative mod- 
els: a simple power-law plus gaussian line model, the warm 
absorber ABSORI model (Done et al. 1992; Zdziarski et al. 
1995), and the ionised disc pexriv model. The latter two 
models are included in the xspec package. The Fe Ka line 
parameters that were found earlier were included in these 
fits but were fixed at their best-fitting values. The results 
of these fits are shown in Table ^. The displayed values of 
£ have been corected for the different energy range in their 
definition, and so can be compared with the values in Ta- 
ble []. The ABSORI models were run with the temperature of 
the absorber fixed at 10 6 K, and with solar Fe abundance 
(adjusted to match the abundances used in the RF models). 
Similarly, the pexriv models had the disc temperature fixed 
at 10 6 K, the inclination angle fixed at 30 degrees, and solar 
abundances (after correction). The power- law cutoff energy 
for the pexriv models was fixed at 100 keV, in agreement 
with our models. 



4.2 Interpretation of the spectral fitting results 

Studying the results listed in Tables gives rise to the 
following comments: 

(i) The RF ionised disc model provides an adequate fit to 
the ASCA data in five out of the six sources examined. The 
only source that the model does not fit is NGC 4051, which 
contains a warm absorber and is well known for its spectral 



complexity (Guainazzi et al. 1996). An attempt to fully ac- 
count for all the features in the spectrum of NGC 4051 is 
beyond the scope of this paper. 

(ii) In general, the alternative models presented in Table [| 
have difficulty fitting the data with realisitic values of the 
parameters. These models seem to be able to fit the data 
only with a large value of the photon index, F. If F is fixed 
at the value determined from the 3-10 keV fits, then the 
fits are quite poor and often results in an extremely small 
Fe Ka EW. This is because of the unrealistically strong Fe K 
edge one finds in the ABSORI and pexriv models (e.g, Ross 
et al. 1999). The RF models determine a value of F that is 
comparable with the ones determined from the high energy 
fit, and so does a good job matching the intrinsic continuum 
shape. 

(iii) The RF results show that the accretion discs in these 
sources might be relatively highly ionised with £ > 1500; 
however, the values determined by ABSORI and pexriv are 
usually too low and poorly constrained. Furthermore, these 
models can only reach a maximum £ of 3200, and so are not 
capable of fitting highly ionised sources. 

(iv) The RF fit for the high luminosity source PKS 0558- 
504 is not well constrained in R — £ space (see the lower- 
right panel in Fig. This is because it is difficult to tell 
the difference between a highly reflective, highly ionised disc 
and a moderately ionised disc with a low reflection fraction. 
In Table |^ we list fits for this source when the reflection 
fraction is free and when it is fixed at unity. Both cases 
give good fits to the data. We will use the parameters from 
the fit when R was free in later discussions, but the above 
degeneracy should be kept in mind. One possible way to 
break this degeneracy is to search for ionised edges and lines 
in the soft X-ray band. These features should be present 
in a spectrum from a moderately ionised disc, even with 
a low reflection fraction. They would not be present at all 
from a highly ionised reprocessor. Obviously, the data would 
need to have extremely good signal-to-noise to detect these 
features and rule out alternatives such as warm absorbers. 

(v) The four sources which have well constrained RF fits 
all have inclination angles between 15 and 30 degrees which 
contradict some earlier results (Ark 564: Turner, George & 
Netzer 1999). These are consistent with the inclination an- 
gles of broad- line Seyfert 1 galaxies (Nandra et al. 1997a). 
Fits to PKS 0558-504 result in large inclination angles. How- 
ever, the PKS 0558-504 Fe Ka line is quite noisy, so those 
inclinations are quite suspect. 

(vi) None of the Fe Ka lines have centroid energies that 
are consistent with coming from highly ionised iron. How- 
ever, this does not mean that ionised Fe is not present in the 
data, because the lines will be redshifted due to relativistic 
blurring effects. Again, the noisiness of the PKS 0558-504 
Fe Ka line did not allow precise measurements of the line 
parameters. 

(vii) Aside from PKS 0558-504, all the sources are best 
fit with a disc emissivity law that is proportional to r -2 ' 5 , 
consistent with the results of Nandra et al. (1997a). Also, 
most of the sources are consistent with the emission ocurring 
inwards of 10 r g . However, it should be emphasized that we 
didn't fit for a and rw„, we just considered four different 
cases. 
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Table 2. Best fit parameters of the RF ionised disc models. All values are quoted in the rest frame. Y is the photon index, £ is the ionisation 
parameter as defined in Equation [l], R is the reflection fraction, r m i n is the inner radius of the emitting 
exponent for the disc emissivity law, and i is the inclination of the disc to the line of sight in degrees. 
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5 DISCUSSION 



The chief advantage of the RF model over the pexriv model 
is that it fits the Fe Ka line and the continuum simultane- 
ously, as they are both computed together. Since the Fe Ka 
line is often the most significant deviation in an otherwise 
featureless power-law spectrum, this ability to fit both com- 
ponants together is a significant improvement. This is best 
illustrated by noting that the RF models find a value of F 
that is consistent with the value computed using only high 



energy data, whereas models like pexriv and ABSORI find 
significantly larger values. 

We can now use our results to compare with ones de- 
rived from simpler models. One of these results is an appar- 
ent correlation between the reflection fraction R, and the 
photon index F, such that objects with larger photon in- 
dices have larger reflection fractions (Zdziarski, Lubihski & 
Smith 1999). In Figure § we plot R versus T for the best 
fitting blurred RF models (excluding NGC 4051, but in- 
cluding the hee-R fit of PKS 0558-504). This figure shows 
that any correlation between R and F must be very weak. 
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Figure 1. Results from fitting TON S f80 with the RF ionised disc model. The top- left figure shows the best fit 
unblurred model: V = 2.48, log£ = 3.30, and R = 0.922. The top-right panel shows the best fit relativistically 
blurred model when i = 30 degrees: T = 2.46, log£ = 3.58, R = 1.22, r m i n = 6 r g , and a = —2.5. The bottom left 
panel displays the ratio of the data to the best blurred model from the top right panel. Here, only the SO data (no 
symbols) and the G2 data (open symbols) have been rebinned and plotted. Finally, the bottom right panel shows 
the 68, 90 and 99 per cent confidence contours around the best fitting blurred model shown in the top right panel. 
The position of the best fit model is denoted by the plus sign. 



However, it should be noted that the correlation pointed out 
by Zdziarski et al. (1999) was mainly for sources with hard 
spectra; they did not have many sources with V > 2. 

It was mentioned in Section ^ that the constant den- 
sity RF models might not be inconsistent with the results of 
the variable density models of Nayakshin et al. (2000) when 
r > 2. Evidence for this interpretation is presented in Fig- 
ure ^, which plots logf of the best fit blurred RF models 
(when the inclination was fixed at 30 degrees) of our sam- 
ple versus the equivalent width of the Fe Ka line that was 
computed from the high energy power-law plus gaussian fits. 
This figure shows a general trend that more highly ionised 
sources have larger Fe Ka equivalent width. This is a typ- 
ical prediction of constant density models (e.g. Matt et al. 
1993): as the top of the disc becomes more ionised, emission 
from highly ionised Fe at 6.7 and 6.9 keV will become vis- 
ible and this will increase the equivalent width of the line. 
Recall that most of our sources have relatively low 2-10 keV 
luminosities (Table [l]), and so we should not observe (and, 
indeed, we do not) a decrease in Fe Ka EW with luminos- 



ity which is seen in higher luminosity sources (the so-called 
X-ray Baldwin effect; Iwasawa & Taniguchi 1993; Nandra et 
al. 1997b). If our sources all had flat spectra (F ^ 2) then 
the trend seen in Figure |^ would be cause for concern; how- 
ever, only one of our sources has F < 2 (plotted with the 
open symbol in Figure ^) and this source has the smallest 
ionisation parameter in our sample, consistent with the re- 
sults of Nayakshin et al. (2000). The rest of the sources have 
T > 2 and have larger ionisation parameters. Therefore, we 
conclude that our results are not inconsistent with the re- 
sults of Nayakshin et al. (2000) despite our use of constant 
density models. 

One other limitation of the models used in this paper 
is the restriction to solar abundance. There is no reason a 
priori to expect that the abundances in accretion discs to be 
solar, but, then again, there is no strong evidence that they 
are not solar. Changing the abundances has the greatest ef- 
fect on the Fe Ka line, but, given the poor signal-to-noise in 
most of the observations presented above, disentangling this 
effect from the other ones would be quite difficult. Hope- 
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Figure 2. Results from fitting PKS 0558—504 with the RF ionised disc model. The top-left figure shows the best 
fit unblurred model: T = 2.27, log£ = 3.69, and R = 0.165. The top-right panel shows the best fit relativistically 
blurred model when i = 30 degrees: T = 2.27, log£ = 3.68, R = 0.198, r mi „ = 10 r g , and a = -2.0. The bottom left 
panel displays the ratio of the data to the best blurred model from the top right panel. Here, only the SO data (no 
symbols) and the G2 data (open symbols) have been rebinned and plotted. Finally, the bottom right panel shows 
the 68, 90 and 99 per cent confidence contours around the best fitting blurred model shown in the top right panel. 
The position of the best fit model is denoted by the plus sign. 



fully, higher signal-to-noise observations by XMM-Newton 
and Chandra will help constrain the abundances in accre- 
tion discs. 



6 CONCLUSIONS 

We have fit ASCA observations of six NLS1 galaxies with 
the ionised disc models of RF, and found that five of the 
sources are well fit by such a model. The only source that was 
not well fit by the ionised disc model was the complicated 
source NGC 4051. The ionisation parameter and reflection 
fraction of the highest luminosity source, PKS 0558-504, 
were not well constrained by the models. This is because 
the models suffer from a degeneracy between high R, high £ 
discs and low R, low £ discs. Observations with higher signal- 
to-noise in the soft X-ray band might be able to break this 
degeneracy. 

The four well constrained sources all have inclination 
angles typical of broad-line Seyfert 1 galaxies (between 15 



and 30 degrees), and are consistent with most of the emis- 
sion coming from inside 10 r g . We find only a very weak or 
nonexistent correlation between the reflection fraction and 
the photon index. 

While our models make the rather large simpilfying as- 
sumption that the density of the disc is constant, because 
most of the sample sources have steep spectra and will there- 
fore show ionised features in their reflection spectra (pre- 
dicted from the non-constant density models of Nakashin 
et al. 2000), our results are probably not unreasonable. In- 
deed, the source with the flattest spectrum has the lowest 
ionisation parameter in the sample and a small Fe Ka EW. 

These results are the strongest evidence yet that NLS1 
galaxies might posess ionised discs. This improvement comes 
about because the RF models include the Fe Ka line emis- 
sion, can handle large ionisation parameters, and we fit a 
large energy range (typically 1-10 keV). Another ionised 
disc model, pexriv, was usually the best out of the alterna- 
tive models considered, but had difficulty fitting the contin- 
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Figure 3. Results from fitting NGC 4051 with the RF ionised disc model. The top-left figure shows the best fit 
unblurred model: T = 1.82, log£ = 3.39, and R = 2.0. The top-right panel shows the best fit relativistically blurred 
model when i = 30 degrees: T = 1.85, log£ = 3.198, R = 1.21, r m i n = 6 r 9 , and a = —2.5. The bottom left 
panel displays the ratio of the data to the best blurred model from the top right panel. Here, only the SO data (no 
symbols) and the G2 data (open symbols) have been rebinned and plotted. Finally, the bottom right panel shows 
the 68, 90 and 99 per cent confidence contours around the best fitting blurred model shown in the top right panel. 
The position of the best fit model is denoted by the plus sign. 



uum and typically underestimated the ionisation parameter. 
The pexriv model is not appropriate to fit possibly highly 
ionised sources like NLS1 galaxies. These results bolster the 
claim that the unique properties of NLS1 are driven by rel- 
atively low mass black holes accreting at a high fraction of 
their Eddington rate. 
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Figure 5. Results from fitting Mrk 335 with the RF ionised disc model. The top-left figure shows the best fit 
unblurred model: T = i.98, log£ = 3.34, and R = 0.678. The top-right panel shows the best fit relativistically 
blurred model when i = 30 degrees: T = f.98, log£ = 3.25, R = 0.888, r m j n = 6 r g , and a = —2.5. The bottom left 
panel displays the ratio of the data to the best blurred model from the top right panel. Here, only the SO data (no 
symbols) and the G2 data (open symbols) have been rebinned and plotted. Finally, the bottom right panel shows 
the 68, 90 and 99 per cent confidence contours around the best fitting blurred model shown in the top right panel. 
The position of the best fit model is denoted by the plus sign. 
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Figure 6. Results from fitting PG 1244+026 with the RF ionised disc model. The top-left figure shows the best 
fit unblurred model: Y = 2.43, log£ = 3.91, and R = 1.44. The top-right panel shows the best fit relativistically 
blurred model when i = 30 degrees: V = 2.42, log£ = 3.86, R = 1.89, r m i n = 10 r 9 , and a = —2.5. The bottom left 
panel displays the ratio of the data to the best blurred model from the top right panel. Here, only the SO data (no 
symbols) and the G2 data (open symbols) have been rebinned and plotted. Finally, the bottom right panel shows 
the 68, 90 and 99 per cent confidence contours around the best fitting blurred model shown in the top right panel. 
The position of the best fit model is denoted by the plus sign. 
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Table 3. Best fit Fe Ka line parameters of the sample galaxies. All values are quoted in the rest frame. V 
is the photon index, Ep e is the energy centroid of the Gaussian fit in keV, cr is the width of the Gaussian 
also in keV and EW is the equivalent width of the line measured in eV. 
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Figure 7. A plot of R vs. T using the data from the best fitting blurred RF model with the inclination fixed at 30 degrees. If there is a 
correlation it is very weak. 
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Table 4. Fit parameters from the alternative models considered. All values are quoted in the rest frame. T is the photon index and EW 
is the equivalent width of the Fe Ka line in eV (the line energy and width are fixed with the values shown in Table ^). £ is the ionisation 
parameter (in units of erg cm s — 1 ) of the warm absorber (for the ABSORI models) or the ionised disc (for the pexriv models). The value of 
£ has been corrected to the same energy range as the RF models. Njj is the column density of the warm absorber in units of 10 22 cm~ 2 . 
R is the reflection fraction given by the pexriv models. The inclination angle of the disc has been fixed at 30 degrees for these models. 
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Figure 8. Plot of log£ from the best fitting blurred RF models (with the inclination fixed at 30 degrees) versus Fe Ka EW, excluding 
NGC 4051. The source with T < 2.0 (Mrk 335) is plotted with an open symbol. 
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